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We propose the possible detection of broken mirror symmetries in correlated two-dimensional 
materials by elastotransport measnrements. Using linear response theory we calculate the “shear 
conductivity” dehned as the linear change of the longitudinal conductivity axx due to a shear 

strain Cxy. This quantity can only be non-vanishing when in-plane mirror symmetries are broken 
and we discuss how candidate states in the cuprate pseudogap regime (e.g. various loop current 
or charge orders) may exhibit a finite shear conductivity. We also provide a realistic experimental 
protocol for detecting such a response. 


I. INTRODUCTION 

Phases of matter in solids are often empirically distin¬ 
guished by their transport properties. Electrical trans¬ 
port measurements probe long wavelength properties of 
the system, and in metallic phases, they are sensitive 
to electronic excitations near the Fermi level. In addi¬ 
tion, these measurements often exhibit singular features 
at phase transitions, thus indicating the onset of broken 
symmetries. However, aside from a few instances such 
as the anomalous Hall effect in ferromagnets, the exact 
form of the broken symmetry is not usually evident from 
a transport measurement. Experiments that aid in di¬ 
rectly identifying subtle forms of broken symmetry are 
therefore invaluable in the study of strongly correlated 
electron materials. 

Motivated by such considerations, we study the linear 
change of electrical transport coefficients in the presence 
of applied strain which we refer to as a “shear conductiv¬ 
ity” . If such a linear response is present, it is necessarily 
encoded in a fourth-rank tensor. As we discuss below, the 
shear conductivity is a binary indicator of point group 
and mirror symmetry breaking and retains its character 
as a transport coefficient in probing the dynamics of the 
quasiparticles (or lack thereof) near the Fermi level. By 
contrast, ordinary transport coefficients, which are sec¬ 
ond rank tensors, are at best indirect markers of such 
transitions. 

Specifically, we discuss how a linear change in longi¬ 
tudinal conductivity axx due to applied strain Cxy can 
only occur if certain vertical mirror plane symmetries 
are broken. When such symmetries are absent a re¬ 
sponse of this sort is no longer forbidden, and is there¬ 
fore generically finite. While our considerations here are 
based solely on symmetry and are therefore quite gen¬ 
eral, we are primarily motivated by the cuprate super¬ 
conductors, where a variety of broken symmetry phases 
are likely present in the pseudogap regime of hole-doped 
materials. Several candidate order parameter theories 
have been proposed for the pseudogap regime including 


current-loop phasei^^^, d-density wave phase^^, v arious 
forms of charge ordei!^®^, electron nematic phase^^^EIl^ 
and pair density wave state^l^f^^ to name just a few. 

Here, we have focussed on two phases that break 
among other symmetries, point group and mirror sym¬ 
metry. These are the variants of phases with loop current 
order as well as those with charge order. Our key result 
is that there is a finite and measurable shear conductivity 
in both states. Furthermore, we predict a parametrically 
higher shear conductivity response in the orbital current 
loop phase near its onset temperature, when compared to 
the response of charge ordered states. Our analysis was 
inspired by an elegant set of experiment that have 
utilized transport measurements in the presence of strain 
as a probe of nematicity. Our goal is to generalize these 
experimental protocols to help uncover more subtle pat¬ 
terns of symmetry breaking. We also note that while our 
focus here is on electrical transport coefficients, the sym¬ 
metry considerations discussed below apply equally to 
other measurements, such as ultrasound attenuatioiP^, 
which also involves the determination of a fourth rank 
tensor. 

This paper is organized as follows: we first review the 
simple symmetry considerations which lead to a non¬ 
vanishing shear conductivity in Sec. |TTj We then discuss 
how this quantity is actually evaluated in the framework 
of the Kubo formula, before performing explicit calcula¬ 
tions for model charge ordered systems and loop current 
phases in Sec. |HI| Next, in Sec. |IV] we discuss how the 
response is trained in macroscopic crystals, before clos¬ 
ing by discussing realistic experimental protocols for the 
measurement of the shear conductivity in Sec. |Vj 

II. SYMMETRY CONSIDERATIONS 

We define the shear conductivity as the elastoconduc¬ 
tivity tensor component Txx,xy = daxx/d^xy, which de¬ 
scribes the change of the longitudinal DC conductivity in¬ 
duced by a shear strain of the crystal. The tensor Txx,xy 
can be non-vanishing only when each of the symmetries 
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is broken: (i) reflection about the ccz-plane, which has a 
normal vector along y: Uy (ii) reflection about the yz- 
plane, and (iii) the combination a(^x=y) * C' 4 - Here, 
^{x=y) denotes reflection about the vertical {x = —y)z 
plane and C 4 ,, a fourfold rotation about the principal z- 
axis. This follows from the fact that under each of these 
symmetry operations, is even while Cxy is odd. More 
generally, in the presence of any of the 3 symmetries men¬ 
tioned above, vanishes if the indices contain an odd 
number of x or y. An example is T^y^xx-, which represents 
the change in the Hall conductivity due to a longitudinal 
strain. On the other hand, inversion symmetry imposes 
no restrictions Tij^ki- Moreover, Onsager’s reciprocity 
theorem dictates that = T where M 

is odd under time-reversal (the final pair of indices is 
unaffected since strain is a symmetric time-reversal in¬ 
variant tensor). Lastly, we remark that reciprocity does 
not relate and Tn^kj- Further details on symmetry 

properties of Txx,xy and the full elastotransport tensor in 
a tetragonal system are provided in Appendix 


III. ELASTOCONDUCTIVITY IN MODEL 
SYSTEMS 


A. Calculation of elastoconductivity 


Having disposed of generalities, we now compute the 
elastoconductivity tensor from an explicit microscopic 
model, using the Kubo formula, and taking into ac¬ 
count the effects of strain. We consider a non-interacting 
Hamiltonian, which captures the appropriate broken 
symmetry, and takes the form 

"H = ^ ^fc.a/3cI.aCfe,/3 = X! (1) 

k,O 60 k 

where a,/3 can be spin, band or other quantum num¬ 
bers, = (Cfe ijCfc 2 i ■ ■ ■) denote the fermionic creation 
operators, and k is the crystal momentum. Anticipat¬ 
ing the physical systems that we will apply our results to 
in the next section, we restrict ourselves to 2D; however, 
the formalism generalizes in a straightforward way to 3D. 
Without strain, the Kubo formula give^^ 
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where Ak{E) = ^\Gk{E + i/ 2 Tsc) - Gk{E - j/St-sc)] 
is the spectral functiorP^ and f{E) is the Fermi distri¬ 
bution function. As explained in the Appendix the 
application of a strain e leads to a change of the Bra- 
vais lattice vectors {fli} —>■ {(l -I- e) a^}, which can easily 
be implemented in a tight-binding approach. This has 
two main effects on a tight-binding HamiltoniarPSl: (i) 
the tight-binding hopping parameters may change since 


they depend on the distance between the atoms in gen¬ 
eral, and (ii) the momenta are modified according to 
fc —(1 -I- e)fe. As a consequence, the Hamiltonian is 
modified as Hk —>■ where the superscript s in¬ 

dicates the modified tight binding parameters, and the 
Brillouin zone is also altered accordingly. 

After introducing strain in ([^, a coordinate transfor¬ 
mation p = (1 -I- e)fc effectively undoes (ii) while in¬ 
troducing a Jacobian into the expression for a^x- As a 
result, the DC conductivity in the strained crystal can 
be written as 
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where the momentum integration spans the unstrained 
1st Brillouin zone of the lattice. The shear conductiv¬ 
ity ^xx,xy'! c&'ii b© computGci from (jsj) by sotting ^xx — 
Cyy = 0, expanding and extracting the linear coefficient 

^xx ~ '^XX + ^xx.xy^xy + G{e^y). 

We now apply the formalism outlined above to study 
two candidate phases relevant to the pseudogap regime 
of the hole doped cuprates: (i) a particular form of 
charge order that breaks mirror symmetries, and (ii) two- 
dimensional loop current phases proposed first by Varma. 


B. Shear conductivity as a probe of charge order 

Despite the remarkable recent experimental progress 
in identifying (generally short range correlated) charge 
order as a ubiquitous member of the cuprate phase 
diagram,f22H3l] j^any fundamental questions of symme¬ 
try remain. Cand idate st ates such as predominantly d- 
wave charge ordeili^lMlMI and criss-crossed stripeJ^^ do 
break the mirror symmetries required to produce a finite 
shear conductivity, while conventional s-wave checker¬ 
board charge density waves, or unidirectional striped 
would produce no such response. A possible detection of 
broken mirror symmetries is signihcant as it reveals the 
presence of robust, thermodynamic, long range ordered 
phase in the pseudogap regime. 

A model calculation can be done by considering the 
most general Hamiltonian which includes charge and 
bond density waves (CDWs and BDWs) in both x- and y- 
directions. This captures the essential mirror-symmetry 
breaking physics, and more complicated candidate states 
(e.g. d-wave charge order) can be constructed from this 
model. At mean held level, the density wave Hamiltonian 
we consider is given by 

S'HdW = E (I^Ri^RpRi+ E + -f H.c.) 

Ri ^ OL—x,y -* 

( 4 ) 
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FIG. 1. Left: A real space picture of a square lattice in the 
presence of a CDW and BDW in both x and y directions. 
Darker points indicate a higher electron density and darker 
lines an increased hopping between neighboring sites. Right: 
A real space picture of the Cu02 plane showing the 02 planar 
loop order phase which produces a finite shear conductivity. 
Orange circles represent copper sites, green circles are oxy¬ 
gen sites, and arrows indicate the direction of spontaneonsly 
formed microscopic current loops. 


where (j)R. = (j)x cos{Qx ■ Ri) + cos{Qy ■ Ri) describes 

the CDW and = ^x/y cos{Qx/y ■ Ri) the BDW. We 
omit the spin degrees of freedom which will only double 
the overall conductivity, and for simplicity we consider 
period two density waves, i.e. \Qx/y\ = ■ With these 

choices, a candidate striped state has (say) finite (j)x and 
Aj; but (f>y = Ay = 0, whtle a state with both s- and 
d-wave charge order in both x and y directions has finite 
(f>x = <Pyi along with finite = — Ay. With similar gen¬ 
eralizations, any candidate commensurate charge order 
can be included in this model Hamiltonian, and while 
the specific choice of ordering period will affect the fol¬ 
lowing symmetry considerations, we consider the above 
an appropriate minimal model to calculate the shear con¬ 
ductivity. 

Before proceeding with a formal calculation, it is use¬ 
ful to anticipate our results using symmetry arguments. 
From the real space picture in Fig. it is easy to see 
that the vertical mirror plane symmetries ax, and ay are 
broken only when finite period 2 charge and bond DW 
fields in both x and y directions are present. Finally, the 
product ax=y * C 4 . is also a broken symmetry, since the 
state is not invariant under four-fold rotations. Thus, we 
only expect a non-vanishing shear conductance for this 
system if all charge and bond fields are finite. Moreover, 
because a mirror reflection is implemented by the change 
in sign of any of the order parameters above, and Txx,xy 
much change sign under vertical mirror plane reflections, 
we expect that the shear conductivity is in general an odd 
power of each order parameter. The lowest order contri¬ 
bution would therefore be linear in each order parameter, 
i.c. F oc (l)x(l)yAxAy 

To apply this analysis to the cuprate like systems, 
we use a single band effective next-nearest neighbor 
tight-binding dispersion Ck = — 2 t[cos{kx) + cos(fcy)] + 
At' cos{kx) cos{ky) — /r of the Cu square lattice with pa¬ 
rameters t' = 0.45t and /i = —0.65t (where t is the near¬ 
est neighbor hopping for the unstrained system). The 
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FIG. 2. Ratio Vxx.xyj of shear conductivity over longi¬ 
tudinal conductivity for checkerboard charge and bond order 
(scattering time Tsc = 10/t). Blue dashed curve as a function 
of (j)x for a finite and fixed value of the other three order pa¬ 
rameters: (j)y ~ Ax = — Ay = 0.11 and red curve as a function 
of Ax for similarly fixed (f>x = 4>y = — Ay = 0.1 1 . 

shear strain does not affect the nearest-neighbor hop¬ 
ping at 0 [txy), but does modify the next nearest neigh¬ 
bor hopping. On the other hand, both the charge and 
the bond orders are unaffected at this order in strain. 
Rewriting the Hamiltonian in the four band basis of 
the new (four site) unit cell, and expanding Eq. ([^ 
to linear order ^xy in tbe limit T —>■ 0 then leads to 

^ XX — ^ XX T ^xx^xy^xy — Px 2PyCxy with 
"K j^a 

= J (^Retr 

— 7Tj2a 

( 5 ) 

where the momentum integral is over the new reduced 
Brillouin zone. The results of the calculation are shown 
in Fig. where we plot the ratio Txx,xy/crxx as a func¬ 
tion of the charge order (f>x and A^,. Consistent with our 
expectations based on symmetry, the response vanishes 
when any of (l)x,(l>y. Ax, or Ay vanishes, since mirror sym¬ 
metries are partially restored when this happens. Close 
to the onset temperature of such charge order, we expect 
the response to be small, as it is proportional to a high 
power of the order parameter fields. However, there is 
no reason to expect a small response well below the on¬ 
set temperature, where the order parameters strengths 
are generically of order unity, in the case of long-range 
charge order. 

C. Shear conductivity as a probe of loop current 
order 

While it is now clear that charge ordered states ex¬ 
ist within the pseudogap regime in many cuprates, a 
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FIG. 3. Ratio /(yxx. for the loop current states 0i 

(brown dashes) and 02 (blue) as function of the order pa¬ 
rameter magnitudes R/t, where t the tight-binding hopping 
parameter of the copper-oxygen bond. The response is para¬ 
metrically larger than for charge order since it is quadratic in 
the magnitude of the small, dimensionless order parameters, 
R/t, and not quartic as in the case of charge order. 


more vexing issue is the nature of the pseudogap it¬ 
self. The pseudogap onset temperature T* is typically 
higher than the charge ordering temperature, with a pri¬ 
mary experimental signature being the suppression of 
the electronic density of stateP^"^, and the onset of 
Q = 0 magnetisrrpSMU Such magnetic order is consis¬ 
tent with predictions first made by VarmgP^, who argued 
that various loop current states are responsible for the 
pseudogap phenomenology. Recent theoretical worlP^HlIl 
has suggested that orbital current phases could also ex¬ 
plain magieto-optic effects observed in the pseudogap 
regim^^®2!. Here we discuss how planar versions of such 
a phase can be detected via shear conductivity measure¬ 
ments. 

We have considered the so called 0i and 02 planar 
loop current states described in Ref. HI As Fig. [^demon¬ 
strates, the 02 state breaks all the required symmetries 
{^xi ^{x=y)*C 4 ) to produce a shear conductivity, while 
the related 0i state with currents on all four diago¬ 
nal bonds preserves the mirror planes Gx and Gy and so 
should not produce a shear conductivity response. 

For a formal calculation, we work with a three band 
model that includes the Cu d-orbitals as well as the 
oxygen px- and Py-orbitals and consider the mean-field 
Hamiltonian with order parameter R of the loop current 
(formally an anapole or toroidal moment, which is a polar 
vector Following the same procedure as before, we 

implement strain into this model by modifying the tight 
binding parameters and shearing the Brillouin zone. The 
results are shown in Fig. and, as expected, the system 
shows a finite shear conductivity Txx,xy in the symmetry 
broken state 02 , but not in 0i. 

We note that unlike the response to charge order, 
^xx,xy is even in the current loop order parameter, R. 
This follows from the symmetry of the 02 state, where 
the state with order parameter R is related to its time 


reversed partner by a C 2 rotation about the z axis (see 
Fig. 0. Because a C 2 rotation leaves the tensor Txx,xy 
invariant, we have the equality ^xx,xy{R) = ^xx,xy{~R)- 
Therefore, the response must be proportional to an even 
power of R. We stress that the shear conductivity is 
likely to be considerably larger in the Varma current 
loop phase than in the charge ordered states studied in 
the previous section since it is proportional to a smaller 
power of the order parameter R near the onset tempera¬ 
ture, where a Landau-Ginzburg treatment remains valid. 


IV. TRAINING THE SHEAR CONDUCTIVITY 
RESPONSE 


Having demonstrated that Txx,xy is finite when the 
requisite point group symmetries are broken, it is natu¬ 
ral to ask whether the shear conductivity is measurable 
in macroscopic crystals where domains are inevitably 
present. Indeed, for both the charge and loop current 
order parameters, the response from different domains 
will cancel, as we discuss below. Nevertheless, because 
the response involves a product of order parameters (ie. 
it is proportional to a composite order parameter), it is 
in fact possible to train this composite order parameter 
(and hence the shear conductivity) across the putative 
phase transition into a mirror symmetry breaking phase. 

One can anticipate how such training is achieved based 
on general symmetry considerations: in crystals where a 
point group symmetry is present in the symmet¬ 
ric phase (the corresponding space group is PA/mmm), 
^xx,xy behaves like the i? 2 g {dxy) representation of 
and so must be proportional to a composite order param¬ 
eter with this symmetry. We can therefore train domains 
of such a composite order parameter by applying a shear 
strain, Cxy, through a symmetry breaking transition. Be¬ 
low, we discuss the formal symmetry considerations that 
allow such training for the loop current phase (where 
the absence of translation symmetry breaking makes the 
analysis simpler), and then the charge ordered phase. 

First, in the case of the 02 loop current state, there 
are four possible domains: two time reversed partners 
with currents on the x = y diagonals (with order pa¬ 
rameters denoted as ±i2“), and another two with cur¬ 
rents on the opposite diagonal {±R’’). The responses 
from domains with currents on opposite diagonals can¬ 
cel, because they are related by a rotation which sends 
^xx,xy(^R ) t T'yy yxi^R ) = ^xx,xy(^R ); where the 
second equality follows from G(^x=y) mirror symmetry of 
the state (see Appendix A2). Nevertheless, because the 
responses from time reversed partners on a given diag¬ 
onal are equal, a macroscopic response is possible if we 
can bias domains to be of a single orientation. This is 
easily accomplished by cooling through the loop current 
transition in the presence of an externally applied i? 2 g 
strain field (a shear strain e^y), which couples to the ori¬ 
entational component of the loop current order in the 

free energy like SF = Xexy where A is a 
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FIG. 4. Experimental setup for the detection of shear resis¬ 
tance similar to Ref. |20l and m] 


coupling constant. 

Next, in the case of the charge order considered in 
Sec. |IIIB[ there are two possible domains of the compos¬ 
ite order parameter which if both present 

will lead to a vanishing shear conductivity. Neverthe¬ 
less, while these individual order parameters cannot be 
trained, their product preserves translation symmetry 
while only breaking point group symmetries and so can 
once more be trained by application of an appropriate 
strain. We can determine the appropriate symmetry 
channel as follows: the composite order parameter ipxAx 
is translation invariant for the case of Qx = (tt, 0). This 
pattern breaks the ax mirror symmetry, and upon C 4 
rotations, is transformed into the composite order pa¬ 
rameter ipyAy. Thus, 4>xAx and /pyAy have Px and Py 
symmetry, and so transform like components of the two 
dimensional representation of D^h- The product of 
these two order parameters is therefore in the i? 2 g rep¬ 
resentation of and there is an allowed coupling to 
shear strain in the free energy, with a term of the form 
6 F = Xexy(t>x4>yAxAy. Applying shear strain and cooling 
through the charge ordering transition should therefore 
induce a finite shear conductivity response. 


V. SHEAR RESISTANCE MEASUREMENTS 


Having computed Txx,xy for various ordered phases, we 
now discuss the experimental protocols required to mea¬ 
sure such a response. Here, we build on recently used ex¬ 
perimental method4i21121l for determining specific terms 


the elastoresistivity tensoJ^ map,-yS = 




£->■0 


which describes the strain-induced change in the normal¬ 
ized resistivity. The technique, which was originally em¬ 
ployed to reveal the presence of nematic fluctuations, can 
be generalized in the following way. 

A rectangular sample with principal axes a, b oriented 
along its edges is glued with a relative angle 9 onto a 
piezoelectric stack in order to measure the longitudinal 
elastoresistance (Fig. [^. Applying a voltage to the piezo 
leads to a stress and therefore tensile strains txx and 
^yy ~ ~^p^xx (where Vp is the Poisson ratio of the piezo¬ 
electric stack) in the coordinate system of the piezo stack. 
Since the sample is mounted with a relative angle on the 
piezo element, the pure tensile strain of the piezo stack 


translates into a combination of tensile and shear strain 
in the coordinate system of the crystal axes a, b for the 
sample. 

Measuring the change of the resistance {^^/r)xx — 
{^P/p)xx to the applied strain for several angles 9, 
the angularly antisymmetric component of the response 
is given by 





V R 


{Vp + 1 ) 


XX 

m. 
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{-9) 
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xx,zz^ zz,xy 
Czz.zz 


^xx^xy 


( 6 ) 

sin(20) • Cx 


where are components of the elastic stiffness ten¬ 

sor. For tetragonal or orthorhombic symmetry, Czz,xy = 
'm‘xx,xy = 0, and the odd contribution vanishes; con¬ 
versely, a nonzero is proof of broken mirror 

symmetries. An equivalent approach involves using litho¬ 
graphic methods to prepare samples with different axes 
orientations (see Appendix]^. 


VI. DISCUSSION 

We have demonstrated the use of a higher rank tensor 
— the shear conductivity Txx,xy — as a sensitive probe 
of broken point group symmetries. In the context of the 
cuprates, the onset of such a response at the charge order¬ 
ing temperature would unambiguously demonstrate that 
charge order breaks all in-plane mirror reflection symme¬ 
tries; recall these are discrete symmetries which are less 
susceptible to disorder and can survive as true long range 
ordered states. However, it is possible that a current loop 
state onsets at a higher, pseudogap temperature scale, in 
which case the detection of Txx,xy at T* supports pro¬ 
posals of the 02 version of this Q — 0 magnetic state. In 
either scenario, any identihcation of broken mirror sym¬ 
metries is significant as it specifies a true thermodynamic 
phase in the pseudogap regime, and places constraints on 
a pntative quantum critical point which may exist near 
optimal doping in cuprate materials. 

To summarize, we have considered transport coeffi¬ 
cients in the presence of strain that are direct indica¬ 
tors of mirror symmetry breaking. We have studied the 
shear conductivity in the context of two broken symme¬ 
try phases that have been proposed to occur within the 
pseudogap regime of hole-doped cuprates and have dis¬ 
cussed experimental protocols to measure such response. 
Our analysis can be generalized to study other compo¬ 
nents of the strain conductivity tensor Tijki- 
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Appendix A: Symmetry considerations 
1. General properties of tensorial response functions 

Given a Hamiltonian % that is invariant under the point group transformation O : r ^ r', a, physical response 
function Rxi,...,x„, which is a tensor of rank n (e.g. resistivity, elasticity, etc.) is constrained according to 

Rxi ,....Xn - • ■ • ^Xn,yn ^yi ,...yn ' 

This condition typically restricts several elements of R to vanish; conversely, breaking some point group symmetries 
relaxes these constraints, allowing the new responses to serve as evidence of broken symmetry. In the following, we will 
restrict ourselves to the case of (quasi) two-dimensional materials as the various cuprates, iron-based superconductors 
or heavy-fermion compounds mostly possess a tetragonal symmetry in the high-temperature regime. Nevertheless, it 
is straightforward to generalize these considerations to more general crystal classes. 

Consider now a state that is invariant under reflection about say the xy-plane : {x, y, z} —>■ {x, y, —z}). We find 
the condition that 


R 


■Xi 


= i-l)^^R 


, where is the number of z-indices (and likewise for ax 
z-axis leads to the conditions 


and ay symmetry). Fourfold rotation symmetry about the 


R 


X\ ,...,Xn 


{-l)^-Rx, 


x-i^y’’ 


and 


Rxi,...,x„ — { 1) Rxi,...,x.„,\^^y- 

Finally, a mirror symmetry a(^x=y) (about the {x = —y)z plane) implies the condition 

Rxi,...,x„ = Rxi,...,x^\^^y- 

Of particular interest to us is the shear conductivity component Txx,xy = daxxldcxy for 2D systems, which describes 
the change of the longitudinal DC conductivity induced by a shear strain of the crystal. Based on the above symmetry 
considerations, this can be non-vanishing only when reflections about the x- and the y-axes - ay, ax - as well as the 
combination a(^x=y) * where a(^x=y} denotes reflection about the vertical {x = —y)z plane and C 4 a tt/2 rotation 
about the z-axis, are broken. Following these lines, it is straightforward to construct the explicit forms of higher-rank 
tensors as the elastoresistance or the elastic stiffness systems with arbitrary symmetries, see e.g., Section [E} 


2. Properties of for planar loop currents 


We now discuss properties of the shear conductivity within a given loop current state, i.e. a single domain which 
corresponds to one configuration of the loop current order parameter R. The magnetic point group symmetry within 
the loop current state is isomorphic to D 2 i^, with the following group elements in addition to the identity: (72(z)T, 
C 2 {x = y)T, C 2 {x = —y), iT, a^, a(x=-y)T, and a{x=y'), where C' 2 (z) denotes a 180 degree rotation about the z axis 
etc., T denotes time reversal, i is inversion, and a[x=y) denotes mirror reflection as before etc. Note that these are 
the group elements for the configuration shown in Fig. These group operations place the following constraints on 
the shear conductivity in the presence of a c axis directed magnetic field, R^'- 


C2{z)T, iT: 

r XX,xy 

[R,Hz) = 

- ^xx,xy{Rj Rz) 

(A2) 

II 

<b 

1 

r XX,xy 

[R,Hz) = 

~ ^yy,yx{Rj ~Hz) 

(A3) 


r xx,xy 

[R,Hz)-- 

~ ^yy,yx{Ri Hz)- 

(A4) 


while dz places no constraint. For clarity, note that G 2 T and iT independently imply the identity of Eq. |A2[ and 
similarly for the subsequent equations. 



3. Properties of for charge order 


In the case of the charge order described in Sec. |IIIB[ the point group symmetry is much more restricted. In the 
case of equal magnitude charge and bond order parameters as drawn in Fig. with \(j)x\ = \(j)y \ and |Aa;| = |Aj,|, the 
point group symmetry is C 2 v, with the group elements C 2 {x = —y), ai^x=y)i and dz in addition to the identity. These 
enforce the following constraints on the shear conductivity in a c axis directed magnetic field: 

^^2(2: = y')^ d'(^x—y) ■ ^xx^xyi^z) — ^yy,yx{ Hz) (-^^) 

while dz places no constraint. Alternatively, when the magnitudes of the order parameters are different, with 
\(j)y\ and |Aa;| ^ |Ay| the point group symmetry is reduced to Cs, with only two group elements: the identity, and 
d'x=y This places no constrains on transformations of Txx,xy 


Appendix B: Linear Response theory of elastoconductivity 


In this section we derive the general expression for the elastoconductivity using Kubo’s Linear Response theory. 
The presented formalism is suitable for any kind of bilinear Hamiltonian with the form 

Hk.o.pcl ^Ck,p = Yy (Bl) 

k,<y^ k 

with corresponding Matsubara Green’s function Gk(ioJn) = — Hk\ ^■ Following MaharPSIthe longitudinal DC 

conductivity for such a multiband system can be calculated to be 



1st B.Z. 



*^Retr 

at 


(^fc(e) 


dHk. 

dkx ' 


(B2) 


where ifc(e) = ^ [Gk{e + i/2Tsc) - Gk{e - i/2Tsc)] is the spectral weight of the system, e is an energy and /(e) 
is the usual Fermi-Dirac distribution function. In order to yield a finite result of this transport expression, it is 
necessary to introduce a finite lifetime for the fermionic quasiparticles, which can originate from electron-electron, 
electron-phonon or electron-impurity interactions. Here, we just assume Tsc ^ 0, and neglect other processes (e.g., 
vertex corrections). Thus, the scattering rate I/t^c ^ Ep (where Ep is the Fermi energy) can be viewed simply as 
an effective parameter of the theory. 


We now focus on two spatial dimensions and consider a system with a Hamiltonian (Bl|, to be described by a 
tight-binding model 


Ri a,/3 


Ci^/3 + 


E 

<.Ri,Rj '> a,0 


E + h-c.] + E E 

<<Ri,Rj» 


oc(3 t 
ij ^i,OL^3 


+ h.c.] -|- ..., 


(B3) 


where Ri are the positions of the atoms of the lattice (Bravais vectors ai, 02 ) and a describes for example possible 
orbital or spin degrees of freedom. The matrix/hopping elements {G)ff are in general dependent on the distance 
between the atoms, which can be parameterized as 


=[e{\Ri - Rj\)]‘'^ = ^ 


(B4) 


where R^ are the equilibrium positions of the atoms. Without strain, the Fourier space representation of (B3) is just 
given by 


= EE4,aCfc./3 

k a.,j3 


-I- 2 t°^P [cos(ai • k) + cos(a 2 ■ fc)] 
-I- 2(t')“^ [cos([ai -I- 02 ] • k) + cos([ai - 02 ] • fc)] -t-... 


k 


(B5) 
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with the spinor = (cfc,i, Ck, 2 : ■ ■ describing the multiband system. We now apply a general strain to the crystal, 
which has the form 


^XX ^X 

£xy Cyi/y 

The strain changes the Bravais lattice vectors ai, 02 of the system to 

a^j = {i + i)aj, 3 = 1,2 


(B6) 


(B7) 


and therefore also the positions of the atoms —>■ i?® = (1 + e)Ri. The Hamiltonian in the presence of strain is 

therefore given by 


k a,/3 


73 + Y, cos([(i + e)a,] ■ k) 

i=i.2 


i=± 


2_^ ^ k^n+e)k^ k- 


(B8) 


Note that there are now three major changes if we compare (B 8 I with (B5): i) Due to (B7), the Brillouin zone is 
strained, which is indicated by ii) the Hamiltoniair matrix 


'^(l+e)fc ^(i+(.)k 




(B9) 


has modified tight-binding parameters depending on the strain e; and iii) the quasimomentum k —> (1 -|- e)fc takes 
into account the change of the Brillouin zone. Inserting now our strained Hamiltonian (B8) into the linear response 
theory we find 


strained 
1st B.Z. 


drk 


de^^^Retr 


de 




\ 2n 

^^{l+e)k ^ 

dkx 


Tre^ f f 

det(i -be) J (27r)2 J 

1st B.Z. —cxo 


oe 


f)fT‘ 

7l^(e)[(l + e..)^+e,,^ 


(BIO) 


where in the second line we substituted p = (1 -I- e)fc so that the integral over the original (unstrained) Brillouin zone 
is recovered. In the following we will always differentiate between the momentum k defined in the strained Brillouin 
zone and p in the unstrained Brillouin zone (see Fig.j^. 


1. Next-nearest neighbor tight-binding dispersion for square lattice in the presence of shear strain 


In order to make the connection between the generalized tight-binding Hamiltonian in (B51 to the systems considered 


in this paper, we now calculate the dispersion of a simple square lattice in next-nearest neighbor approximation and 
in the presence of a pure shear strain 


e = 


0 


^xy 

0 


(Bll) 


The Bravais vectors of a square lattice are ai = e^, 02 = Gy, 


where we set the lattice constant a = 1. The nearest- 


neighbor hoppings = t (and also the charge and bond order parameters) do not change linearly in e^jj, by a pure 


shear strain since the distance between neighbor atoms is |(1 -b \ = 1 -b 

neighbor hoppings 


1. In contrast, the next-nearest 


[t'{t 


. ± e«)]^ = ^ 


(BI2) 
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unstrained sheared 



FIG. 5. Unstrained and sheared Bravais lattice and corresponding first Brillouin zone of square lattices. 


obtain a linear ea;j,-dependence and we get the sheared dispersion 

Ep = 2t[cos(pa;) + cos(pp)] + 2t'(l - \/2-f'e^y) cos{px + Py) 

+ 2t'(l + V2j'ej:y) cos{p^ -py) - p.. (B13) 

Note that p is, as mentioned in the last section, the momentum in the unstrained Brillouin zone which has to be used 
in formula (BIO). 


Appendix C: Shear conductivity in cuprates with charge order 

1. Effective Hamiltonian for bidirectional charge order 

A bidirectional charge {(f)) and bond (A) density wave of period A S N on a square lattice (which describes our 
quasi two-dimensional material) is described by the commensurate fields 


^Ri = ■ cos{Q^ ■ R,) + (fy ■ cos{Qy ■ R,), 

= ^x/y ■ cosiQ^/y ■ Ri), 


(Cl) 


where <f)x/y,^x/y G R and Qx/y = W^xly Note that symmetric (i.e., “checkerboard”) order with cfx = 4’y and 
Ax = Ay, also preserves the &(^x=-y) symmetry. We consider here the effective square copper lattice (Bravais lattice 
(J-x = Oy = ^y) describing the material’s dispersion (note that we have ignored spin degrees of freedom which 
simply contribute an overall factor of 2 in the conductivity). The effective Hamiltonian reads 


= [t + +1' ^ c]ci-bh.c. 

Ri OL—x,y Rj—Ri-\-aa Rj—Ri-\-ax:tay 


(C2) 


where we use the parameters t = l,t' = 0A5t,p = —0.65t. Transforming this Hamiltonian into Fourier space leads to 


'H- = ^ekclck+Y^ H cos[(fc-f-^) • a„] c^g^Cfc-f h.c.Y 

k (y.—x,y ^ ^ ' 


(C3) 


with the free dispersion 


Efc = — 2t[cos(fca;) -I- cos{ky)] + 4t' cos(A:a;) cos{ky) — p. (C4) 

Due to the density waves the periodicity of the crystal changes from a = 1 to A, so we have downfolded our original 
Brillouin zone 


kx/y e [-tTjTt] -)► kx/y G [-7r/7V,7r/A], 


(C5) 
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FIG. 6. Fermi surface of the original (left) and down-folded (right) Brillouin zone for charge and bond order {N = 2) with 
gaps opening at the hot spots of the Fermi surface that are connected by = Tre^/y. Note that the figure on the right is 

the reduced Brillouin zone shown here as the upper right corner of the original Brillouin zone. 


SO that we end up with an effective band Hamiltonian 


k 


(C6) 


where '$i. = (c^, Cfc+g ,.. .)^. From this point on we will restrict ourselves to the case of = 2, but the following 
arguments and steps are valid for arbitrary N. For N = 2 we have the spinor 'Fj. = {ck,Ck+Q^,Ck+Qy,Ck+Q^+Q^)^ 
and the Hamiltonian 


Hk = 


( Cfc 

Xx(fc) 

Xl{k) 

V 0 


Xx{k) 

efc+Qx 

0 

X*y{k + Qx) 


Xy{k) 

0 

^k+Qy 

X%{k + Qy) 


0 \ 

Xy{k + Qx) 
Xx{k+Qy) ’ 

efc+Q.+Qj, / 


(C7) 


with 


Xx/yik) — 4^x/y “b ‘^^xjy^ COs[(/2 -f Qyj^ ■ Q-gj/yj — ^xjy ^^^xjy ^^^\kx/y “b 

In Fig. we show the original cuprate Fermi surface and the down-folded Fermi surface for the period 2 charge and 
bond density waves. 


2. Shear conductivity 


We now discuss how to apply the shear strain (Bll) on the density wave system. The change of the Bravais-lattice 
“r/y = (1 + will also lead to a change of the reciprocal lattice vectors x cby/x) finally also 


the ordering vector Q® = T 6® ^y. This means that 


Qa ' ®/3 — Qc 


2tt 

iV’ 


(C9) 

the nearest neighbor hoppings are 
by a pure shear strain and the same argument should therefore hold for the charge and 
bond order (jjx/y, ^x/y since their ordering vectors are aligned with the x and y axis. This means that the strained 


remains invariant under the application of the strain. As we saw in Section B 1 
not affected linearly in 


Hamiltonian of (C7) is given by 


h; = 


X*xiP) 

X*y{p) 

0 


■Xp) 


0 


Xy{p) 

0 

P+Qy 


X*y{P+ Qx) X*x{.P+ Qy) 


Xy{P + Qx' 
Xx{P + Qy 


(CIO) 
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with the strained dispersion e® defined in (B13). We now rewrite the formula for the elastoconductivity in the case 
of pure shear strain 


g fmol i 7re"^ 

'^xx — ~ 1 _ 2 

^ ^xy 


1st B.Z. 


k/2 

-77/2 


(fp 


de^^Retr 

oe 


{^pi^)Jp,x 


(Cll) 


where the integral is performed over the reduced Brillouin zone and we have defined the current vertex as 


P = 

^ p,x 


dm 


dm 


dpj; dpy ■ 

Expanding the spectral weight and the current operator to linear order in the shear strain e^y we have 

Alie) 


(C12) 




P = P 

^p,x ^p,x 






^xy — Ap{e) + dy4p(e) 


py^J ^xy 




f)f 


(C13) 


^xy — Jp^x H” ^Jp^x ‘ ^xy- 




Inserting this in (Cll I and expanding again, we find the DC conductivity 

{^Ap{^Jpx) 


(7'P'r' — TT e 


/ ^ 


/-7r/2 (277) 


de 


(C14) 


and the shear conductivity 


F — 

XX,xy — ^/i D 


r' d^P f , dfie)^ 


d-Ap{^c)Jp xAp{c)Jp x ~\~ Api^c^S Jp^x-Ap(^6^ Jp^; 


(C15) 


A numerical investigation shows that Txx,xy is insensitive to 7 ', so we henceforth set 'y' = 0. In this case we can 
finally simplify 


^XX - Qxi 


^xx^xy — ‘^Qyj 


(C16) 


with 


Q ct — 


LXJ 

' d'^P f , 5/(e)„ 

' deyAdRetr 

oe 


2 

->■ 776 


-7r/2 (277)y 

d2p 


1-1^12 


( 277 ) 


Re tr 


1 \ 

dHr, 


^Px ^Pi 


(C17) 


where in the end we considered the zero temperature limit. This formula is then numerically evaluated to find the 
shear conductivity numerically, as discussed in the main text. 


Appendix D: Shear conductivity in cuprates for current loop order 

Following Ref. [2] we calculate the shear conductivity for the Varma current loop states 0/ and 0//. From the 
minimal model of the copper oxides that takes into account the dx^-y^ orbitals of the copper atoms and the Px,Py 
orbitals of the oxygen atoms we can write down the mean-field Hamiltonian of the time-reversal breaking fields as 
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(spin degrees of freedom are introduced later) 


y-011 


k 


k 


'dk \ 

Cd-^i 

[2tpd - Re^* 

Px,k I 

h.c. 



' \ h.c. 

h.c. 

'dk \ 

Vcd - M 

‘2‘tpdSx{k) 

Px,k j 

h.c. 

-t 

,Py,k) 

' V h.c. 

h.c 


AtppSxik^Syih^ 


AtppSx{}^^Sy{k) 

-M 



(Dl) 


where tpd and tpp are the next nearest neighbor hopping parameters of the copper-oxygen and oxygen-oxygen bonds, 
td the on-site energy of the copper orbitals, /x the chemical potential, R ■ the current loop order parameter of the 
Varma theory (with (p = ±7r/2) and we defined Sx/y{k) = sin(fc • ax/y), Cx/y{k) = cos(fc • ax/y) with ax/y = a ■ exjy 
The creation/annihilation operators are defined in such a way that the current vertex is again just given by 


k 


dj^ 

dk 


(D2) 


Introducing the shear as in the previous section and assuming a linear shear-dependence of = tpp{l — jexy) 

and t^i^xy) = tpp{l + jexy) (where arro ws schematically indicate the direction of oxygen-oxygen bonds on the 
diagonal), we can expand expression (BIO) linear in strain. For the numerical investigation we use the parameters 
tpp = 0.3tpd,ed = —0.7tpd,/x = 1.5tpd,T = 10/tpd,^ = l.Stpd- As discussed in the main text, the 0/ state, which does 
not break the and ay mirror symmetry, does not produce a finite shear conductivity — Txx,xy = 0. In contrast, 
the 0/7 state breaks the requisite mirror symmetries and hence shows a finite shear conductivity Txx,xy ^ 0 which 
increases with the order parameter R. Note that it is the 0// state (and several symmetry related partner states 
involving out of plane current loops to apical oxygens) which is both a stronger theoretical and experimental candidate 
for a pseudogap state. 


Appendix E: Experimental Setups 


1. Explicit forms of tensors for broken tetragonal symmetry 


Let us now consider (quasi) two-dimensional materials with an intrinsic Oz symmetry, e.g., the layered cuprate or 
iron pnictide compounds. If all tetragonal point group symmetries dx, dj,, iT(a;=_y), (74 are broken, a generic rank-4 
tensor only has components with an even number of z-indices, see Sec.|^ In the usual Voigt notation (1,2,3,4, 5, 6) = 
(xcc, yy, zz, yz, zx, xy) we can then write a rank 4 tensor T that is symmetric in the first and last two indices (such as 
the elastoresistance or elastic stiffness)!^ as 


Symmetry: d^ 


/Tn 

Ti 2 

Ti 3 

0 

0 

T^ieX 

T21 

T22 

T23 

0 

0 

T26 

T31 

T32 

T33 

0 

0 

T36 

0 

0 

0 

T44 

T45 

0 

0 

0 

0 

T54 

T55 

0 

\Tei 

762 

Tea 

0 

0 

tJ 


(El) 


Let us now consider what happens if we restore the different point group symmetries d^, d^, d(x=-j,), (74 of the 
tetragonal group. Restoring either dx or dy will lead to the vanishing of all tensor elements with an odd number of 
X-indices (and therefore also for the elements with an odd number of y-indices): 


Symmetries: az,{ax or dy) 


/Til 

Ti2 

Ti 3 

0 

0 


T21 

T22 

T23 

0 

0 

0 

T31 

T32 

T33 

0 

0 

0 

0 

0 

0 

T44 

0 

0 

0 

0 

0 

0 

T55 

0 

V 0 

0 

0 

0 

0 

Tee/ 


(E2) 
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Restoring the a(^x=-y) symmetry relates many 


Symmetries: ^z,^{x=-y) 


Finally, having a and C 4 symmetric system 


Symmetries: az,Ci 


the tensor elements 

in (El I so that 




/Tn 

Ti 2 

Ti3 

0 

0 

Ti 6 \ 



Ti 2 

Tn 

Ti3 

0 

0 

Ti 6 

T — 

^31 

T 31 

T 33 

0 

0 

T 36 



0 

0 

0 

'I'44 'J 45 

0 



0 

0 

0 

T 45 T 44 

0 



VTei 

Tei 

T 36 

0 

0 

Tee/ 

would end up 

with 







/Tn 

Ti 2 

Ti3 

0 

0 

Tie \ 


Ti 2 

Tn 

Ti3 

0 

0 

-Tie 

f = 

T 31 

T 31 

T 33 

0 

0 

0 

0 

0 

0 

T 44 

T 45 

0 


0 

0 

0 

—T 45 

T 44 

0 


\ni 

-Tgj 

0 

0 

0 

Tee / 


(E3) 


(E4) 


Combining the tensor representations (E2)-(E4) one can determine the tensor for an arbitrary combination of the 
considered point group symmetries. E.g. restoring the tetragonal symmetry would lead to the tensor 


Tetragonal symmetry: 


/Tn 

Ti2 

Ti3 

0 

0 


Ti 2 

Tn 

Ti3 

0 

0 

0 

T 31 

T 31 

T 33 

0 

0 

0 

0 

0 

0 

T 44 

0 

0 

0 

0 

0 

0 

T 44 

0 

V 0 

0 

0 

0 

0 

Tee/ 


(E5) 


2. Setup 1 

In this section we describe the two proposed experimental setups to detect tetragonal symmetry breaking in elas- 
toresistance measurements in detail. Earlier experiment^ glued single crystals of Ba(Fe;^_,j.Coa;) 2 As 2 to the surface 
of a piezo stack as shown in Fig. FtI Varying the strain of the piezo crystal exx,£yy = —Vptxx (with Vp the Poisson 
ratio of the piezo stack) Kuo et alT^ measured select admixtures of elastoresistivity coefficients.. Both the resistance 
and the applied strain are tensors of second rank and can be written in the reducible Voigt notation by six component 
arrays 


^ i^xx: ^yy 7 ^zz: ^yz 7 ^zx 7 ^xy ) 


^ {^xx7 ^yy 7 ^ 


yy 7 ^zz 7 ^yz 7 


\T 

^zx ? ^xy) 


(E 6 ) 


Defining the relative change of the resistance due to the applied strain as Sri = we can write down the relation 


Sri = rh^jej, 


(E7) 


with the elastoresistance tensor m. Note that although m is a tensor of rank-4 (as it relates two rank-2 tensors), 
the Voigt notation simplifies this tensor to a matrix form. Consider gluing a rectangular sample, which is cut along 
the a, b crystal axes, onto the stack with a relative angle of 0 as shown in Fig. The sample is described by the 
elastoresistance tensor rh and the elastic stiffness C (both tensors are in the basis that is shown along the a, b crystal 
axes). Assuming that we have a system w hich breaks all of the tetragonal symmetries, we can write down the elastic 
stiffness and elastoresistance according to (El) a^ 


/Cn 

C12 

Ci 3 

0 

0 

Cie\ 


/mil 

mi 2 

mi 3 
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0 

mie\ 

^12 

C'22 

C 23 

0 

0 

C2e 


W 21 

TO22 

m23 

0 

0 

m26 

Ci 3 

<^23 

C'33 

0 

0 

C' 3 e 


W 31 

TO32 

W 33 

0 

0 

TO36 

0 

0 

0 

C44 

C45 

0 

^ Tn — 

0 

0 

0 

m44 

m45 

0 

0 

0 

0 

C45 

C 55 

0 


0 

0 

0 

m54 

TO 55 

0 

VCie 

C2e 

C 36 

0 

0 

Cee/ 


Vmei 

TOe2 

me3 

0 

0 

Wee/ 
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FIG. 7. The first example of an experimental setup, where samples which have the same alignment of their crystalline a — b 
axes are glued at several different angles to a piezoelectric stack. 


Upon restoring the tetragonal symmetry, this simplifies to 


C = 


1^11 

C/12 

C/13 

U 

u 

u 



mil 

777 1 2 

777 1 3 

U 

u 


Cl2 

Cii 

^13 

0 

0 

0 



77112 

77711 

777 1 3 

0 

0 

0 

Ci 3 

^13 

^33 

0 

0 

0 



77131 

77731 

77733 

0 

0 

0 

0 

0 

0 

0*44 

0 

0 


, m = 

0 

0 

0 

77744 

0 

0 

0 

0 

0 

0 

^44 

0 



0 

0 

0 

0 

77744 

0 


\ 0 0 0 


\ 0 0 0 0 


(E9) 


Comparing (E 8 ) with (E9) we see that, as stated earlier, the breaking of the in-plane mirror symmetries allow for a 
broad range of new responses (e. g., t he elastoresistivity coefhcient mie = mxx,xy ^ 0 ). 


We now consider the relation (E71 for the rotated sample and in the absence of in-plane mirror symmetries. We 


align our coordinate system axes with the rectangular crystal edges and therefore along the a, b crystal axes. In this 
system, the elastoresistance and elastic stiffness are given by the expressions in (E 8 | and the strain Ccrystai seen by 
the crystal is related to the strain in the coordinate system of the piezo by a rotation 

^crystal — RcXgR ^piezo • (ElO) 

Here, we defined the relation between the real and the engineering strain R = diag(l, 1,1,2, 2, 2) and the rotation 
matrix for the 6 -component arrays similar to Ref. [20j : 
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ae = 


cos^( 0 ) 
sin^ { 0 ) 
0 
0 
0 


cos^{0) 

0 

0 

0 


0 0 0 

0 0 0 

1 0 0 

0 cos( 0 ) — sin( 0 ) 
0 sin( 0 ) cos{6) 


y— cos( 0 ) sin( 0 ) cos( 0 ) sin( 0 ) 0 0 


0 


2 cos( 0 ) sin( 0 ) ^ 
—2 cos( 0 ) sin( 0 ) 

0 

0 

0 

cos^( 0 ) — sin^( 0 ) y 


(Ell) 


The strain in the piezo system is given by Cpiezo = —Vp, —Vzi0),O,O,O)£xx- Note, that due to the crystal axis 
rotation the Poisson ratio ^z(0) of the crystal in z direction now depends on the rotation angle 9 (see Fig. [^. The 
angular dependence can be derived by considering the stress-strain relation Tcrystai = Cccrystai and the condition 
Ts,crystal = Tzz,crysted = 0 since there is no applied stress in z direction of the crystal. One finds that 

Ci 3 (cos^(0) — izsin^(0)) -I- C 23 (sin^(0) — jzcos^(0)) — + 1) sin(20) 




C- 


33 


The resistance change is therefore given by 

Sr = TOCcrystai = mRaeR~^ < 

such that the measured aiai-component is given by 


-piezo? 


(E12) 


(E13) 


Sri{9) = (^) = ^mii [cos^( 6 >) — z/pSin^( 6 *)] -I-mi 2 [sin^(0) — VpCOS^(9)\ 


- miQ{vp + 1 ) sin(26l) - 77 x 131 /^( 6 *) 


(E14) 


Looking at the odd contribution we finally find 


AR 


odd 

{9) = 


AR 


(9)- 


AR 


(-9) 


— i^p + 1 ) 


^xx,zz^zz,xy 


Cz 


^xx,xy 


sin(26*) 
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3. Setup 2 

The second proposed setup is similar to the one described in the previous section, but here we start from a thin 
film crystal and use lithography to cut out the rectangular sections as shown in Fig. The advantage of this method 
is that aligning the crystal by gluing it on the piezo element with high precision is a rather difficult task, whereas in 
using lithography, one can cut out the samples very precisely. In contrast to the other experiment, the crystal axes 




FIG. 8. A second possible experimental setup, where crystals with different alignment of their crystalline axes can onto the 
piezoelectric stack. 


a, b are aligned with the principal axes of the piezo stack. Therefore, the Poisson ratio in z direction does not depend 
on the angle of the cut. If we consider the coordinate system to be along the a,b (and thus also x^y) axis, what is 
measured is effectively the rotated resistance 


dr{9) = agSro = agm ■ Cpiezo- 
The measured longitudinal resistance is therefore given by 


(E15) 


6ri{B) = (^) = ^cos^(^) [wii - vmi 2 - + sin^(6') [m2i - i'm 22 - ^^2^23] 

- sin(26») [mgi - vmQ 2 - wea] ^ ^xx , 


(E16) 


and the antisymmetric part 

odd 

(6») = - sin(2d) [mgi - vmQ 2 - i^zWea] ^-xx, 

XX 


/Ai? 

[ji 


is again proportional to sin(20) and only present for broken mirror symmetries. 


(E17) 



















